% 指定 Excel 文件路径和名称
filename = 't-s-v电源直连.xlsx';
sheet = '25.158';
% 指定要读取的数据范围，例如：'A1:B10'
dataRange = 'A223:C1988';
% 使用 xlsread 函数读取数据
[data, headers] = xlsread(filename, sheet, dataRange);
% 将数据存储到 MATLAB 的变量中
throttle = data(:, 1); % 第一列数据
motor_speed = data(:, 2); % 第二列数据
voltage = data(:, 3); % 第三列数据
[p,S] = polyfit(throttle,motor_speed,6);
[fit_ms,fit_delta] = polyval(p,unique(throttle),S);

Fs = 100; % 采样率（每秒采样点数）
T = 1/Fs; % 采样间隔
% 创建采样数据
L = length(data(:, 2)); % 数据长度
t = (0:L-1)*T; % 时间向量
% 计算频谱
Y = fft(motor_speed); % 对信号进行FFT变换
P2 = abs(Y/L); % 双边频谱
amp1 = P2(1:L/2+1); % 单边频谱
amp1(2:end-1) = 2*amp1(2:end-1); % 去除负频谱并加倍幅值
% 创建频率向量
f1 = Fs*(0:(L/2))/L;

figure(1);
plot(throttle,motor_speed,'o',unique(throttle),fit_ms,'r-',unique(throttle),fit_delta,'b-');

figure(2);
plot(f1,amp1);
